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■ Abstract 

. QCD equations for generating functionals are solved at coinciding 

I momenta of particles. As a result, the relations for q-particle corre- 

^ ■ lation functions at equal momenta are obtained. They are directly 

OO . connected with the previously derived results about the factorial and 

I cumulant moments of multiplicity distributions. It is predicted that 
the correlations at coinciding points decrease, first, as a function of 

. the rank q, and then start oscillating when represented in a form of the 

I ratio of cumulant to factorial correlations. In particular, the cumulant 

' function of the fifth rank should become negative at the origin. Some 

p • experimental data are briefiy discussed. 

^ : 1 INTRODUCTION 



Multiparticle correlations in hadropro duct ion at high energies are very com- 



■ pUcated both for experimental studies and for theoretical description. Beside 

^ rather sophisticated dynamics, the problems due to the huge number of avail- 

able variables enter the game. The choice of any particular variable leads to 
some specific correlations being emphasized (see, e.g., Ref. |]l|). 

The multiparicle correlation functions integrated over momenta give rise 
to the moments of multiplicity distributions. Some special features of the mo- 
ments have been revealed from the solution of QCD equations for generating 
functions (for the review, see Ref. 0). 

Surely, the dependence of the correlation functions on particle momenta 
provides more detailed information. Unfortunately, their study in the 3q- 
dimensional phase space is not an easy task. That is why up to now we have 
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got experimental data on 2-particle (and, sometimes, on 3-particle) correla- 
tion functions only plotted as functions of a particular variable (usually, it is 
the relative (pseudo) rapidity or 4- momentum squared). 

Theoretical approach asks for the equations for generating functionals. 
In QCD they were written rather long time ago P, However, no general 
solution of them has been found. It is common to rewrite those equations in 
the form of cascade equations a'la Altarelli-Parisi, and to find their solution 
for inclusive distributions or for some specific correlations P, H, ^, |[. 

We shall show how to get some information about correlations directly 
from the equations for generating functionals using the experience we got 
from the solution of corresponding equations for the generating functions [0]. 
We demonstrate it with the illustrative example of the relations between the 
correlation functions of various ranks at coinciding momenta. In particular, 
some predictions for the correlations of higher order (yet unavailable in ex- 
periment) are obtained. The more sophisticated case of unequal momenta 
will be considered in a separate publication. 



2 DEFINITIONS AND NOTATIONS 

The main ingredient of any correlation analysis is the g-particle inclusive 
density 

p/-)(k„...,k,)^l (1) 
^ a dui, . . . , duq 

where the invariant momentum phase space is given by 

dcoq = r^-^, (2) 
' (27r)32Eg' ^ ^ 

and the set of indices fi denotes the internal quantum numbers of all q 
particles. We shall omit it in what follows. 

The generating functional of inclusive distributions is defined by the for- 
mula 

°° 1 /■ f 
G[z{k)] =Y,— dui... dUgPgziki) . . . z(kg), (3) 

q=0 1- 

where the integration is over available 3g-dimensional phase space. 
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From d^) one gets 
The normalized factorial correlation function is defined by 

i^.(ki,...,k,)= (5) 

pi(ki)...pi(kj 

The logarithm of the functional is the generating functional of the 
so-called cumulant correlation functions 

"^^'^"•••''^^^ = fe(k,)...^z(k,) '-°- 
The corresponding normalized function is 

In quantum field theory, the functions pq are determined by the whole 
set of Feynman graphs while the functions Kg are given by the connected 
graphs only. In correlation analysis, correspondingly, pg take into account all 
correlations, and Kg correspond to the "genuine" g-particle correlations. 

Integration of these functions over momenta gives rise to factorial and 
cumulant moments of multiplicity distributions. Their generating function 
is described by the same formula as (^) at constant values of the functions 
2;(k) = z = const. 



3 Equations for generating functionals and 
their solution 

The QCD Lagrangian defines the dynamics of hadropro duct ion, and, conse- 
quently, the correlations of particles produced. It determines the relationship 
of generating functionals in the form of non-linear integro-differential equa- 
tions which were proposed a long time ago P, §]. To simplify the equations, 
we shall consider here the gluodynamics where quark degrees of freedom are 



3 



neglected. It is reasonable for qualitative conclusions we are aimed at as 
our previous experience of solution of equations for generating functions 
shows. 

Thus the equation for the generating functional of gluodynamics is 



dG{z{k),Y) 
dY 



dxK{x)^^\G{z, Y+lnx)G{z, Y+\n{l-x))-G{z, Y)]. (8) 

Jo 



Here Y = ln{p6/pQ),p is the initial momentum, 6 is the angular width of the 
gluon jet considered, po =const , 



TT 

0:5 is the coupling constant, and the kernel of the equation is 

K{x) = --{l-x)\2-x{l-x)]. (10) 

X 

It is the non-linear integro-differential equation with shifted arguments in 
the non-linear part which take into account the conservation laws, and with 
the initial condition 

G(z,Y = 0) = l + z{k = p), (11) 

and the normalization 

G{z = 0,Y) = l. (12) 
In particular, at constant z one gets the generating function 

G(z,y) = ^p„(r)(i + ^)", (13) 

n=0 

where P„ is the probability of ?7,-particle events. 

The condition (0) normalizes the total probability to 1, and the condition 
(|TT]) declares that there is a single particle at the very initial stage. The 
corresponding equation for the generating function (|13]) looks quite similar 
and is obtained from eq. (||) at z = const. Its solution was considered in [0]. 
In a similar way, one can get the solution of (§) at the special point where all 
the momenta of the q registered particles (actually, gluons in our treatment) 
are equal i.e. ki = k2 = . . . = kg. Therefore, the momenta differences are 
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chosen equal to zero. It means that we are considering the correlations at 
coinciding points where Ak = 0. Since the correlation functions depend on 
the differences of momenta and on the overall momentum (i.e. on Y) the 
only dependence left in our case is the dependence on Y. 

After Taylor series expansion and differentiation in eq. (^ similar to the 
procedure in [0], one gets the same differential equation 

(InG(F))// = 7o'[G'(F) - 1 - 2/iiG/(F) + h^GffiY)], (14) 

where h = 11/24; /ig = (67 - 6n^)/36 ^ 0.216. 

In general, one can define the anomalous dimension 'y{uJi,as{Y)) by the 
formula which generalizes ones commonly used for average multiplicity or 
single inclusive distributions (see [0, ||) 

fY 

Pg{uj^,Y) = pg(Yo)exp[q dyj{uji, asiy))]- (15) 

JYq 



Substituting ([T5| ) in (|1^) one gets after differentiation for terms with equal 
powers of z^Ui) 



Kq{Ui) Kq{uJi 

nq{Ui) — — 



7^[1 - 2hiq-f{uJi) + h2{q^j^{uJi) + g7'(^i)] 



(16) 



Even though we write down the explicit dependence on uji, it was assumed 
in ([l^ when differentiating over Y that correlation functions are constant in 
Ui i.e. the correlation lengths are much larger than the interval of momenta 
considered. That is why we omit uji in what follows considering it equal to 
zero. Apart from it, the condition qk <C E, where E is the total energy, 
should be fulfilled to ensure no additional dependence on Y. Then all calcu- 
lations are just the same as for the moments of multiplicity distributions in 
Ref.0, and from the normalization condition Hi{0) = 1 one gets the same 
relation between 7(0) and 70 

7(0) ^ 70 - Ih^'o + li^h, - h\)^l + 0(7o^). (17) 

The value of 70 is measured in experiment (it is equal 0.48 at the mass of Z^- 
boson). Therefore it is possible to estimate the values of correlation functions 
at coinciding points according to (0) and (|1^. 
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The main conclusion is that the ratio of cumulant to factorial correlation 
functions at the origin should be approximately equal to the ratio of their 
integrals 

ff,,0).^4£l.^.^,. (18) 



According to results of [0, 0], it predicts that this ratio should decrease fastly 
at g = 2, 3, 4 and reach the negative value at g = 5. At ever higher values of 
q it should oscillate (for the review see 0). 



4 Theoretical results and experimental data 

Our main theoretical statement relates two sets of experimentally accessible 
quantities, namely, moments of multiplicity distributions and the values of 
correlation functions at the origin. It could be verified without any theoretical 
formulae. However, one should clearly recognize possible shortcomings which 
could be of various origin. First, the prediction is done for gluons while in 
experiment we deal with pions mostly. However, the success in predicting 
the corresponding qualitative features of multiplicity distributions (see f^) 
is encouraging. More important problems could be related to hadronization 
effects just at the origin because of resonances, Bose-Einstein correlations etc. 
Moreover, we have now the experimental data only about two- and three- 
particle correlation functions (with rather low precision) while the data about 
moments of multiplicity distributions are much more precise and admit good 
predictions. It appeals to experimentalists. 

At the moment I was able to get the necessary ratios from the data about 
densities Pg provided by NA22 Collaboration using the formulae 

H (0) = ^ = /^3(0)-3pi(0)p2(0) + 2p;(0) 

' P3(0) P3(0) ■ ^ ^ 

Using the data for non-singlediffractive Tr+p events at 250 GeV/c, one gets 

H^^{0) = 0.34; H^^%0) = 0.13; (21) 
H,- -(0) = 0.28; H, (0) = 0.1; (22) 
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H++{0) = 0.25; H+++{0) = 0.07 (23) 

for charged (c), negatives ( ) and positives (+). These values agree quah- 

tatively with those shown in for moments of multiphcity distributions in 
pp data (e.g., UA5 Collaboration data for charged- charged correlations are 
H2 = 0.26; = 0.07) and show the same tendency as the ratios of moments. 
Hardly, more elaborated tests can be done at present stage but one could try 
to check it in Monte Carlo models as well. The most impressive prediction 
of the negative value of the fifth order correlation function at the origin can 
be also verified in Monte Carlo models. 

5 Conclusions 

The solution of the equations for generating functionals in gluodynamics gives 
a hint to the equality of the ratios of integral cumulant and factorial moments 
to the corresponding ratios of the values of cumulant and factorial correlation 
functions at the origin. In particular, it follows that latest ratio should also 
decrease from g = 1 to g = 5 becoming negative, and oscillating at ever higher 
values of ranks q. It would be desirable to verify the prediction using the 
data of same experiments. The comparison to Monte Carlo models is welcome 
too. From the theoretical side it would indicate once again the importance 
of a parameter 57 in QCD (at the inclusive level, it is closely related to 
the exclusive factor nlag widely discussed in vector boson production at 
threshold and in connection with Parke- Taylor amplitudes). Besides, the 
negative values of the ratio at asymptotically high energies would show the 
decline from the negative binomial distribution which predicts the positive 
values of Hg at any q. 
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